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Abstract 



Large-time behavior of solutions to the inflow problem of full compressible 
Navier-Stokes equations is investigated on the half line R+ = (0, +oo). The wave 
structure which contains four waves: the transonic(or degenerate) boundary layer 
solution, 1-rarefaction wave, viscous 2-contact wave and 3-rarefaction wave to the 
inflow problem is described and the asymptotic stability of the superposition of the 
above four wave patterns to the inflow problem of full compressible Navier-Stokes 
equations is proven under some smallness conditions. The proof is given by the ele- 
mentary energy analysis based on the underlying wave structure. The main points 
in the proof are the degeneracies of the transonic boundary layer solution and the 
■ wave interactions in the superposition wave. 
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o5 1 Introduction 

In this paper, we consider an initial-boundary-value problem for full compressible Navier- 
Stokes equations in Eulerian coordinates on the half line R + = (0, +oo) 

pt + (pu)x = 0, 

(pu) t + (pu 2 +p) x = (pu x ) x , x > 0, t > 0, (1.1) 

[p (e + \u 2 )] t + [pu (e + \u 2 ) +pu] x = (k6 x + puu x ) x 

where p(t,x) > 0, u(t,x), 8(t,x) > 0, p(t, x) > and e(t,x) > represent the mass 
density, the velocity, the absolute temperature, the pressure, and the specific internal 
energy of the gas respectively and p > is the coefficient of viscosity, n > is the 
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coefficient of heat conduction. Here we assume that both p and k are positive constants. 
Let v = -(> 0) and s denote the specific volume and the entropy of the gas, respectively. 
Then by the second law of thermodynamics, we have for the ideal polytropic gas 

p = R v ~ 1 9 = Av- 1 exp f ^-l s ] , e(v, 9) = -^—9, (1.2) 

\ R J 7-1 

where 7 > 1 denotes the adiabatic exponent of gas, and A and R are positive constants. 
We consider the initial-boundary-value problem (II. ip with the initial values 

(p,u, 9)(0,x) = (p ,u ,9 )(x) -> (p + ,u + ,9 + ) as x -> +00, inf (p , 0o)(z) > (1.3) 

where p + > 0, u + and + > are given constants. 

As pointed out by [T5], the boundary conditions to the half space problem (11.11) can 
be proposed as one of the following three cases: 

Case I. outflow problem (negative velocity on the boundary): 

u(t,x)\ x=o = u-<0, e(t,x)\ x= o = 6-. (1.4)i 

Case II. impermeable wall problem (zero velocity on the boundary): 

u(t,x)\ x=o = 0, 9(t,x)\ x=0 = 9_. (1.4) 2 

Case III. inflow problem (positive velocity on the boundary): 

u(t,x)\ x=0 = u- > 0, p(t,x)\ x=0 = p-, 9(t,x)\ x=0 = 9_. (1.4)3 

Here all the p_ > 0, w_ and 9- > in (1.4) are prescribed constants and of course 
we assume that the initial values ( 11. 3ft and the boundary conditions (1.4) satisfy the 
compatibility condition at the origin. Notice that in Cases I and II, the density p_ on the 
boundary {x = 0} could not be given, but in Case III, p_ must be imposed due to the 
well-posedness theory of the hyperbolic equation (11.1 1) x . 

In the present paper, we are concerned with the large-time behavior of the solutions to 
the inflow problem (Case III) of the full compressible Navier-Stokes equations (1.1), (1-3) 
and (1.4)3. The large-time behavior of the solutions to the compressible Navier-Stokes 
equations ( 11. ip is closely related to the corresponding Euler system 

{Pt + (pu) x = 0, 
(pu) t +(pu 2 +p) x = 0, (1.5) 
[p(e+i)] t +[pu(e + i)+pu] x = 0. 

The Euler system (11.51) is a typical example of the hyperbolic conservation laws. It is 
well-known that the main feature of the solutions to the hyperbolic conservation laws 
is the formation of the shock wave no matter how smooth the initial values are. The 
Euler system (II. 5p contains three basic wave patterns, that is, two nonlinear waves, called 
shock wave and rarefaction wave and one linear wave called contact discontinuity in the 
solutions to the Riemann problem. The above three dilation invariant wave solutions and 
their linear superpositions in the increasing order of characteristic speed, i.e., Riemann 
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solutions, govern both local and large-time behavior of solutions to the Euler system 
and so govern the large-time behavior of the solutions to the compressible Navier-Stokes 
equations (11. ip . 

There have been a large amount of literature on the large-time behavior of solutions 
to the Cauchy problem of the compressible fluid system (II. ip towards the viscous version 
of the basic wave patterns. We refer to PQ, 0, [5], [7], [S], [II], [13], H3, US], [20], [23], 
[23] and some references therein. All these works show that the large-time behavior of 
the solutions to the Cauchy problem is basically governed by the Riemann solutions to 
its corresponding hyperbolic system. 

Recently, the initial-boundary value problem of fll.ip attracts increasing interest be- 
cause it has more physical meanings and of course produces some new mathematical 
difficulties due to the boundary effect. Not only basic wave patterns but also a new wave, 
which is called boundary layer solution (BL-solution for brevity) [15], may appear in the 
IVBP case. Matsumura [15] proposes a criterion on the question when the BL-solution 
forms to the isentropic Navier-Stokes equations, where the entropy of the gas is assumed 
to be constant and the equation (11.ip 3 for the energy conservation is neglected. The ar- 
gument in [15] for the isentropic Navier-Stokes equations can also be applied to the full 
Navier-Stokes equations f ll.ip . see [3] for details. Consider the Riemann problem to the 
Euler equations (11. 5p . where the initial right state of the Riemann data is given by the 
far field state (p + ,u + ,9 + ) in (II .3p . and the left end state (p_,u_,6L) is given by the all 
possible states which are consistent with the boundary condition (1.4) at {x = 0}. Note 
that to the outflow problem, p_ can not be prescribed and is free on the boundary. On 
one hand, when the left end state is uniquely determined so that the value at the bound- 
ary {x = 0} of the solution to the Riemann problem is consistent with the boundary 
condition, we expect that no BL-solution occurs. On the other hand, if the value of the 
solution to the Riemann problem on the boundary is not consistent with the boundary 
condition for any admissible left end state, we expect a BL-solution which compensates 
the gap comes up. Such BL-solution could be constructed by the stationary solution to 
Navier-Stokes equations. The existence and stability of the BL-solution (to the inflow or 
outflow problems, to the isentropic or full Navier-Stokes equations) are studied extensively 
by many authors, see [3], 0], [6], [10], [15] [18], [21], [25], etc. 

Now we review some recent works on the large-time behavior of the solutions to the 
inflow problem of the full Naiver-Stokes equation ( 11. ip . ( 11.31) . (1.4) 3 by Huang-Li-Shi [3] 
and Qin-Wang pi] . In [21] . we rigorously prove the existence (or non-existence) of BL- 
solution to the inflow problem ( 11. ip . ( 11.31) . (1.4) 3 when the right end state (p + ,u + ,9 + ) 
belongs to the subsonic, transonic and supersonic regions respectively. When (p±, u±,9±) 
both belong to the subsonic region, the BL-solution is expected and the stability of this 
BL-solution and its superposition with the 3-rarefaction wave is proved under some small- 
ness assumptions in |3j. The stability of the superposition of the subsonic BL-solution, the 
viscous 2-contact wave and 3-rarefaction wave is shown in [21] under the condition that 
the amplitude of BL-solution and the contact wave is small enough but the amplitude of 
the rarefaction wave is not necessarily small. The stability of the single viscous contact 
wave is also obtained in [21] if the contact wave is weak enough. It should be remarked 
that the subsonic BL-solution decays exponentially with respect to £ = x — cr_i, which is 
good enough to get the desired estimates. When the boundary value (p_, m_, 9 J) belongs 
to the supersonic region, there is no BL-solution. Thus the large-time behavior of the 
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solution is expected to be same as that of the Cauchy problem and the stability of the 
3-rarefaction waves is also given in [3]. 

In the present paper, we are interested in the stability of wave patterns to the inflow 
problem (ll.ip . (11. 3p and (1.4) 3 when (p_, «_,#_) belongs to the transonic region. In 
this case, a new wave structure which contains four waves: the transonic (or degenerate) 
BL-solution, 1-rarefaction wave, viscous 2-contact wave and 3-rarefaction wave, occurs. 
Due to the fact that the first characteristic speed on the boundary is coincident with the 
speed of the moving boundary in the transonic BL-solution case, the nonlinear waves in 
the first characteristic field may appear, which is quite different from the the regime that 
(p_, m_, 9 J) belongs to the subsonic region in our previous result [21 J, where the waves in 
the first characteristic field must be absent. Here we just assume that the 1-rarefaction 
wave appear in the first characteristic field. Correspondingly, some new mathematical 
difficulties occur due to the degeneracy of the transonic BL-solution and its interactions 
with other wave patterns in the superposition wave. In particular, the transonic boundary 
layer solution is attached with 1-rarefaction wave for all time, so the interaction of these 
two waves should be carefully treated in the stability analysis. 

Because the system (II. ip we consider is in one dimension of the space variable x, it is 
convenient to use the following Lagrangian coordinate transformation: 



p(r,y) dy- pu(r,y) dr 




Thus the system (II. ip can be transformed into the following moving boundary problem 
of Navier-Stokes equations in the Lagrangian coordinates [IB] : 



0. 



(v,u,9)(0,x) = (v ,u ,9 )(x) -> (v + ,u+,0+). 

(v, U, 9){t, X = cr_t) = (f_, U_, #_), U- > 



t > 0, x > a_t, 



as x — > +oo, 



1.6) 



where <r_ := — ^— < is the speed of the moving boundary. 

In order to fix the moving boundary x = cr_t, we introduce a new variable £ = x — a_t. 
Then we have the half-space problem 

f v t — — = 0, 

(A* + y) t - °- + + ^ =k (^), + » ' w 

(v,u,6)(t = 0,0 = (v ,woA) (0 -> (v + ,u + ,9 + ) as £ -> +oo, 
{(v,u,6)(t,£ = 0) = (u_,u_,0_), m_ >0. 

Given the right end state (v + ,u + ,8 + ), we can define the following wave curves in the 

phase space (v, u, 9) with v > and 9 > 0. 

• Transonic (or degenerate) boundary layer curve: 



BL(v + ,u+,e+) := \ [v,u,. 



11 11 

- = -o- = —, (u, 9) e E(u+, 9+] 

V V- 
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where (v+,u+,9+) G Tf rans = {(u, 9)\u = y/Ry9 > } is the transonic region defined 
in (12. 4 p with positive gas velocity and T,(u + ,9 + ) is the trajectory at the point (u + ,9 + ) 
defined in Case II of Lemma 2.1 below. 

• Contact wave curve: 

CD(v + ,u+,9 + ) := {(v,u,9)\u = u + ,p = p + ,v ^ v + }, (1.9) 

• i— Rarefaction wave curve (i = 1,3): 

Xi < X i+ , u = u + - J \i(r), s + ) dr), s(v, 9) = s + | , (1.10) 



Ri(v + ,u + ,9 + ) := < (v,u,9) 



where s + = s(v + , 9 + ) and Aj = Aj(t>, s) is the i— th characteristic speed given in (12. 2p . 
Our main stability result is, roughly speaking, as follows: 

• Assume that (t>_, «_,#_) G BL-R 1 -CD-R 3 (t> + , u + , 9 + ), that is, there exist the 
unique medium states (t>*, m*, 9*) G r^. ans , {v mi u m , 9 m ) and (v *, u*, 9*), such that {v _, m_, 9. 
BL(v*, u*,9*), (v*,u*,9*) G i?i(u m ,w m ,6l m ), (v m ,u m ,9 m ) G CD(v*,u*,9*) and (v*,u*,9*) G 
Ra(v + , m + , 6 I + ), then the superposition of the four wave patterns: the transonic (or degen- 
erate) BL-solution, 1-rarefaction wave, 2-viscous contact wave and 3-rarefaction wave is 
time- asymptotically stable provided that the wave strength 5 = \(v + — f_, u + — u_, 9 + —9J)\ 
is suitably small and the conditions in Theorem 2.1 hold. 

This paper is organized as follows. In Section 2, after giving some preliminaries on 
boundary layer solution, viscous 2-contact wave, rarefaction waves and their superposition, 
we state our main result. In Section 3, first the wave interaction estimations are shown, 
then the desired energy estimates are performed and finally our main result is proven. 

Notations. Throughout this paper, several positive generic constants are denoted by c, C 
without confusion, and C(-) stands for some generic constant(s) depending only on the 
quantity listed in the parenthesis. For function spaces, L P (R + ), 1 < p < oo, denotes the 
usual Lebesgue space on R + . W k ' p (R + ) denotes the k th order Sobolev space, and if p — 2, 
we note H k (R + ) := W k ' 2 (R + ), || ■ || := || ■ || L 2(r + ), and || • || fc := || • || H fc (R+ ) for simplicity. 
The domain R+ will be often abbreviated without confusion. 



2 Preliminaries and Main Result 

It is well known that the hyperbolic system (I1.5P has three characteristic speeds 



AxM) = -^, A 2 = 0, X 3 (v,9) = y^-. (2.1) 

V V 

The first and the third characteristic field is genuinely nonlinear, which may have nonlinear 
waves, shock wave and rarefaction wave, while the second characteristic field is linearly 
degenerate, where contact discontinuity may occur. 

Let 

c(v, s) := y/-v* Pv (v, s) = ^9 =: c(v, 9), M(v, u, 9) := -}^- (2.2) 

c(v,9) 
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be the sound speed and the Mach number at the state (v,u, 9). Correspondingly, set 



c + := c{v + ,9+) = v^T^, M + := M(v + ,u + ,9+) = ^ (2.3) 

c + 

be the sound speed and the Mach number at the far field {x = +00}. We divide the 
phase space {(v, u, 9)\ v > 0, 9 > 0} into three parts: 

tt sub := {{v,u,6) I M < 1 } , 
Ttrans :={(v,u,6) \ M = 1} , (2.4) 

^super ■= {(V,U,9) \ M > 1} . 

Call them subsonic, transonic and supersonic region, respectively. Obviously, if we add 
the alternative condition u > or u < 0, then we have six regions £lf ub , T^; aris , and £lf uper - 

2.1 Boundary layer solution 

When (f_, u-, 6-) G ft+ b U r t + rans , we have 

Ai(u_, 6L) = - ^ Rl6 - < _?h = a _ < o, (2.5) 
hence a stationary solution (V 6 , f/ b , © b )(£) to the inflow problem (11.71) is expected 



( a_V b -U b = 0, 



- + 1 (^) 2 ) £+ (w) 5 (2 - 6) 



l(nf/ ft ,e 6 )(o) = (v 6 ,f/ b ,e fe )(+oo) = (v + ,u+,9 + ) 



where P 6 := p(V b , O 6 ) = We call this stationary solution (V b , U b , b )(O the bound- 
ary layer solution (simply, BL-solution) to the inflow problem (jl.7p . 

From the fact that V b (£) > and w_ > 0, then 



U b u+ u- 

u+>0, — = -t = — = - a _. 2.7 

V° V-L V- 



Thus ( 12.61) is equivalent to (12.71) and the following ODE system 



A 



(u b )' = - a T v b {u b -u + ) + f (e b - e ^v b 

(e b )' = -^v*(e> - e + ) + E±v b (u b - u+ ) + ^v b (u b - u+ )\ (2.8) 
(u b , e b ) (o) = o 6 ) (+00) = 

where p+ := p(f+, 6L|_). 

We can compute that the 

Now we state the existence results of the BL-solution to (12. 8ft while its proof has been 
shown in I2T1. 
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Lemma 2.1 (Existence of BL-solution) |21j Suppose that v± > 0, tt_ > 0, 9± > 

and let 5 b := \(u+ — U-,9 + — 6J)\. If u + < 0, then there is no solution to If 
u + > 0, then there exists a suitably small constant Sq > such that if < 5 b < So, then 
the existence and non-existence of solutions to (12.81) is divided into three cases according 
to the location of (u + ,8 + ) : 

Case I : (u + ,9 + ) G Qf up - Then there is no solution to Ii2.8\) . 

Case II : (u + ,9 + ) G Tf rans . Then (u + ,Q+) is a saddle-knot point to (12 .8j) . Precisely, 
there exists a unique trajectory £ tangent to the straight line 



fMU+(u - u+) - k(7 - 1)(0 - 9+) = 



(2.9) 



at the point (u + , 9 + ). For each («_, 6 J) G 6 + ), there exists a unique solution (U b , b ) 

satisfying 



U\ > 0, 



and 



Jn 



u. 



-,e b -9 + ) 



0(1) 



e*>o, 



rn+l 
°b 



n = 0,1,2,.... 



(2.10) 



Case III : (u + ,9 + ) G ^ ub - Then (u + ,9 + ) is a saddle point to ( 12. 8p . PPrecisely, there 
exists a center-stable manifold M. tangent to the line 



(1 + a 2 c 2 u + )(U B - u+) - a 2 (& B - 9+) = 

on the opposite directions at the point (u + ,9 + ). Here c 2 is one of the solutions to the 
equation 

'Mil - 1 



y 2 + 



M%R-y /<7-l)/ y M$RrfK 







and a 2 = ~ ^x^x' 2 ) w ^ > ^, < ^ are ^ e ^ wo eigenvalues of the linearized matrix 
of ODE (12. 8p . Only when («_,#_) G jM(tt+,#+), does there exist a unique solution 
(U b ,Q b ) cM(u + ,9 + ) satisfying 



d 



— (U b -u + ,e b -9 + ) 



0(l)8 b e-^, n = 0,1,2,.... 



(2.11) 



2.2 Viscous Contact Wave 

If (v-,U-,9-) G CD(v+, u+, 9 + ), then the following Riemann problem 



v t - u x = 0, 

Ut + Px — 0, t > 0, x G R, 

(^0 + §« 2 ) t + (H* = o, 

(t>_, tt_, 6L), x < 0, 



(u,«,0)(O,z) 



(2.12) 
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admits a contact discontinuity solution 



(v, u, 0)(t, x) 



(«_,«_, 0_), 



x < 0, t > 0, 



u+,6+), x > 0, i > 0. 



From [7], the viscous version of the above contact discontinuity, called viscous contact 
wave (V d , U d , d ) (t, x) can be defined by 



ne sim (t,x) 



V d (t,x) 
TTdff T \ — „. , (7-i)«e™(t,x) 



Q d (t,x) = e 



Sim / x 



(2.13) 



+ i? /i- 



(7-l)« 



05 



where 6 sim f ' 
equation 



l+t 



is the unique self-similar solution to the following nonlinear diffusion 



(y-l)np+ /e 



< - ~~ fF 7 

e(t,±oo) = e±. 



le/i» 



(2.14) 



Note that £ = x — cr_£, we have the following Lemma: 

Lemma 2.2. [7J T/ie viscous contact wave (V d ,U d ,Q d )(t,x), (x = £ + cr_t) defined in 
(I2.13P satisfies 

i) ^"(e rf - g ± ) = 0(1)^(1 + t)~f exp (- g^g^i ) , n = 0,1,2, 

ii) t^(t,0 = 0(1)5,(1 + *)- 1 exp(-^g^ 



iii) (^ d ,C/ d ,e d )(t,e = O)-(7;_,u_,0_) = 0(l)5 d e- rf . 

where 5d = \9 + — 9-\ is the amplitude of the viscous contact wave and Cd,c > are 
constants. 

Then the viscous contact wave (V d , U d , d ) defined in (I2.13P satisfies the system 



V d - <r_V£ - Uf = 0, 



R 



(-y> 



(2.15) 



+ H° 



where P d : = p(V d , Q d ) and 



# d = 0(1)^(1 + ^ exp 



c d (z + a-ty 

l + t 



due to Lemma 2.2. 
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2.3 Rarefaction waves 



It is well known that if (t>_, «_,#_) G Ri(v + ,u + ,9 + ), {i = 1,3), then there exist a 
i— rarefaction wave (v r % u ri , 9 n )(x/t) which is the global weak solution to the following 
Riemann problem 



v t - u x = 0, 

Ut + Px = 0, t > 0, x G R, 

(t>_, u-,9-), x < 0, 
0+), a; > 0. 



(2.16) 



Consider the following Burgers equation 
Wt + ww x = 0, t > 0, x G R, 

tuo(x) := w(0, x) = 




x < 0, 

^e - ^ dy, x > 0. 



(2.17) 



+oo 



Here q > 14 is a constant to be determined, and C q is a constant such that C q / y 9 e y d?/ 

Jo 

1. If < then the solution to the above Burgers equation can be expressed by 

w(t,x) = w (x (t,x)), x = x (t, x) + w (x (t,x))t. (2-18) 

Moreover, we have 

• w(t,x) = u>_, if x < W-t. 

• For any positive constant a > and for x > 



= |tu (zo(*, 


X)) 


- w+\ 






r+oo 


= C q (w + ~ 


W- 


) / y q e~ v dy 






J X()(t,x) 






r+oo 


= Cg(w + - 


u;_ 


) / y q e~ 






J x—wo(xo(t,x))t 






r+oo 


< C q (w + - 




) / V q e~ y dy 






J x—w+t 


< C q (w + - 




)e- ao \ if x> 



■o+w+)t (2.19) 

Note that the estimation in (12. 19ft play an important role in the wave interaction estimates, 
which is motivated by [12] and [16J . 

Now the i— rarefaction wave (V Tl , U r \ Q Ti )(t,x) (i = 1, 3) to the inflow problem (jl.7p 
can be defined by 



' \i(V n ,Q ri )(t,x) = w{l + t,x + aJ), 

8(y*,e*)(t,x) = s + = s{v + ,e+), 

fV T i (t,x) 

U n (t,x) = u+- / \i{r),s + )dr). 
9 



(2.20) 



Then the i— rarefaction wave (V r \ U r \ ri )(t, x), (i = 1,3) defined in (I2.20p satisfies the 
system 



Vp -a.Vp -Up = 0, 
Up - a_Up + Pp = 0, 



— <7_ 



T Q n + \{U r *) 2 ^ + (P ri U ri )t = 0, 



(2.21) 



(tm£r*,e*)(U = 0) = («_,« -,9.), 
(y r *,U r *,e r *)(t,0 -)> (v+,u + ,9+) as e^+oo 



where P Ti := p(V ri , Q ri ) . 



Lemma 2.3 i— rarefaction wave (V ri ,U Ti ,O ri )(t,l;), (i = 1,3) defined in (12.201) satis- 
fies 

i) up(t,o>o, (\vp\,\ep\)<cup ; 

ii) For any p (1 < p < oo), there exists a constant C pq such that 

\\(Vp,Up,ei%t)\\ LP < C p min{5 ri ,^(l + t)- 1+1 ^}, 
\\(Vp, Dg, 9- )(t)|| LP < C p min {<y P( , 0^+^(1 + t)- 1+1 /,} ; 

iii) ForVa > ; if £ > [-a. + X 1 (v + ,9+) +2a Q ] (l+t), then &^{{V r \ U r \ 9 ri )(*, €)- 
(v + ,u + ,9 + )}\ < C8 n e-"°\ 71 = 0,1,2,.-. ; 

iv) For £ < [-<t_ + A 3 («-, 0-)] (1 + f), (^ r3 , f/ rs , 6 ra ) - (v_, w_, 0_) = 0; 

v) lim sup \(V ri ,U ri ,G ri )(t,£) - (v ri ,u n ,9 ri )(^-:)\ = 0. 

Remark: The statement iii) is a direct consequence of the (12.191) . 



2.4 Superposition of transonic BL-solution, 1-rarefaction wave, 
2-viscous contact wave and 3-rarefaction wave 



In this subsection, we consider the case that (t>_,u_,0_) G BL-Ri-CD-R 3 (v + ,u + ,9 + ), 
that is, there exist uniquely three medium states (v*, u*, 9*) G Tf rans1 (v m ,u m ,9 m ) and 
(v*,u*,9*) such that (u*,u*,0*) G BL(v-,u-,0J), (u*,u*,0*) G Ri(v m ,u m ,9 m ), (v m ,u m ,9 m ) G 
CD(v*, u*,9*) and (u*, u*,9*) G -Rs(f+, «+,#+)• In fact, three medium states (t>*, «*, 0*) G 
r^ans; ( v m,u m ,9 m ) and (t>*,-u*,0*) can be expressed explicitly and uniquely by the fol- 
lowing nine equations 











< 




Urn 


= w , 


* 

U 


= u + + 



0* 

* ' 

V 



(2.22) 



i?7^^+ tT"^ dry, v*^ L 9* = vT~ L 6+. 



7+1 



*7— 1 n* 



,7-1, 
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Define the superposition wave (V, U, Q){t, £) by 



u 
V / 



(t,0 



V b + V ri + V d + V rs 
U b + U ri +U d + U r ' 3 

Qb + @n + Qd + QT 3 



\ 



^ V* + V m + V* » 



it* + u m + u* 

9* + dm + 



(2.23) 



where (V b , ?7 6 , 6 )(£) is the transonic BL-solution defined in Case II of Lemma 2.1 with 
the right state (v+, u + , 9 + ) replaced by (u*, w*, 0*), (V n , £/ ri , 6 ri )(t, £) is the 1-rarefaction 
wave defined in (12.201) with the states (u_, w_, 6L) and (t> + , w + , # + ) replaced by (f*, «*, 0*) 
and (v m ,u m ,9 m ) respectively, (V^ d , U d , Q d )(t, £) is the viscous contact wave defined in 
(I2.13P with the states (t>_, w_, 6L) and (t> + , w + , 9 + ) replaced by (v m , u m , 9 m ) and (v*,u*, 9*), 
respectively, and (V T3 , U r3 , Q r:i )(t, £) is the 3-rarefaction wave defined in (12.201) with the 
left state replaced by (v*,u*,9*). 

Now we state the main result of the paper as follows. 



Theorem 2.1 (Stability of superposition of four waves) Assume that (t>_, «_,#_) £ 
BL-Ri-CD-R 3 (v + ,u + ,9 + ). Let (V,U,Q)(t,£) be the superposition of the transonic BL- 
solution, 1-rarefaction wave, viscous 2-contact wave and 3-rarefaction wave defined in 
A2.23\) . Then there exists a small positive constant So such that if the initial values and 
the wave strength 8 = \(v + — — u_,9 + — 6L)| satisfy 

8 + \\(v - V Q ,u - U ,9 - 00)11! < 8 . (2.24) 

the inflow problem (1.7) has a unique global-in-time solution (v,u,0)(t,£) satisfying 



(v-V,u-U,9- Q)(t, eC([0,oo); H\R+)), 
(v-V)t(t,OeL 2 (0,oo-L 2 (R+)), 
( U -U,9- 0) c (t, G L 2 (0, oo; H\R+)) . 



(2.25) 



Furthermore, 



lim sup \(v -V,u-U,9 -Q)(t,C>\ = 0. 



(2.26) 



Remark. In Theorem 2.1, we assume that 8 = \(v + — v~,u+ — U-,9 + — 9S)\ is suitably 
small. This assumption is equivalent to the one that the amplitudes of the four waves are 
all suitably small. In fact, from the relations in (12.221) and the facts C/| > 0, U^ 1 > 0, 
U^ 3 > 0, we have 

— V-\ + |0* — 6L| = 0(1) (u* — U-), 
\v m -v*\ + \9 rn - 6*\ = 0(l)(u m - u*), (2.27) 
\v+ - v*\ + \9+ - 9*\ = 0(l)(u + - u*). 

Thus 5b = 0(l)(u* — U-),8 n = 0(l)(u m — u*), 5 r3 = 0(1) (u+ — u*). Due to u m = u* 
by the contact discontinuity curve, we have if 5 is small, then 8b, 8 ri and 8 r:j are all small. 
Furthermore, we have 8 a = \9* — 9 m \ < 8b + 8 ri + 8 r3 + 8 is small. 
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3 Stability Analysis 



3.1 Wave interaction estimates 

Recalling the definition of the superposition wave (V, U, Q)(t, £) defined in (I2.23p . we have 
( V t - a.V^ -Uz = 0, 



U t -a-Ut + Pz = n(!f) + G, 



^ (e t - a_e € ) + put = k(^^ + ^ + h, 



(3.1) 



t (V, U, G)(t, £ = 0) = (?;_, 0_) + (V d , U d , Q d ) (t, £ = 0) - (v m , u m , 6 m ). 
where P := p(V, 6) and 



G= {P - P b - P r > - P d - P r *)t- ^ ~$ - ^ 



f r 3] 



< # = (PUz - P b Ul - P^Ul 1 - P d Uf - P r3 Ul 

(u e r Of) 2 (ui)' 



K ' V yF J + P I 



V 



■w 



(3.2) 



To control the interaction terms coming from different wave patterns, we give the 
following lemma which will be critical in the energy estimate in Subsection 3.3. 

Lemma 3.1 (Wave interaction estimates) 

\V^(V n -v*)\ + \vp{V b -v*)\d£ = 0(l)5 1/s (l + t)- 13/w , 
\V b (V d -v m )\ + \V d (V b - v m ) | d£ = 0(1)5(1 + t)~\ 
\Vl(V r3 -v*)\ + \Vp (V b - v,) | <% = 0(1)5 1/8 (1 + t)- 7/8 , 
\V d (V r ^ -v m )\ + \vp (V d - v m ) | di = 0(l)5e~ c \ 
\V d (V rs -v*)\ + \Vp(V d -v*)\d£= 0{l)5e- ct , 
|yn (yrs - v *)\ + \vp (V'i - v m ) | <% = 0(l)5e" rf , 



R 4 



R.4 



R 4 



R4 



R 4 



(3.3) 



\VlV{\ d£ = 0(1)5(1 + ty 



[ \V*vp\d£ = 0(1)5(1 + t)-\ 
[ \V ( d Vp\d£ = 0(l)5e- ct , 

Proof. First we prove f)3.3P i . that is 



/ \V b Vp\d£ = 0(1)5(1 + t)-\ 
[ \V d Vp\d£ = 0(l)5e- c \ 
[ \VT 1 Vp\dZ = 0(i)6e- et , 



(3.4) 
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• Interaction of transonic boundary layer solution and 1-rarefaction wave: 
Since Vp > and Vg > 0, we have V ri — v * > and — V b > 0. Thus we have 



2 



](l+t) H-oo 



2 " " + / )VP{y,-V h )d£ 

[Jo 7[_ ff _ + M-^)] (1+t) J 5 V ; 

J1 + J2. ' (3.5) 



Note that 

—a- + - 



2 

u- Xi(v m ,9 m ) _ \i(v m ,9m) 
V- 2 t>* 2 

— ^ 1 2 ~~ —M{v*,V*) H 2 

[Ai(t; m ,g w )-Ai(^,g,)]- Al( ^ ,M 



Al(t>m,6m) 

Now we can compute that 



r- g _ + *l(WnO l(l +t ) 







- -_+ > i(^ 8 ") ](i + «) 5 



0(1)11^(^)11^ / 







0(1)<%(1 + t)-8ln(l + 5 ft t) 

0(1)4(1 + *)"^, (3.6) 



and 



00 



= 0(l)6 b (v m - V ri (t,C))\^ = [ a _ , ^jWai j^ 

= O(l)5 b e- CT0 *. (3.7) 

due to the statement iii) in Lemma 2.3 by taking o"o = — Al ^™' gm ) > 0. So the combination 
of ([31]) and ([32]) g ive s O r 
Then we prove (I3.3IW : 

• Interaction of transonic boundary layer solution and viscous 2-contact wave: 
/ \V{(V d -v m ) I + \V*(V b -v*)\d£ 

' c_ t "V 

£l\\Vl(V d -v m )\ + \V*(V b -v*)\dZ 

J 3 + J 4 - 
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We calculate 



T t 

J 3 = [ \Vs b (V d -v m )\ + \V*{V b -v*)\dt 
Jo 

a_t 



0(1)5, j Q exp^ YT^)^ 

0(l)6 d e- ct . (3.1 



Also, we have 



/+oo 
\V*{V d -v m )\ + \V d (V b - u„) | 

• — J 4. J 4. ■ 



We can estimate 



and 



/oo 
a _ t \Vl(V d -v m )\d£ 

= 0(1«(1 + S b tr 2 jH, exp (- Cdi \ + + a ; t)2 ) di 

/oo 2 
exp (-C dV 2 ) drj 
-00 

= 0(1)5,(1 + t)~ 3 / 2 , (3.9) 

/oo 
„JV?(V b -v*)\d£ 

= 0(1)5A(1 + S b t)-\1 + t)~W j~_ t exp (- Cdi \ + + a ; t)2 ) d£ 

= 0(l)5 d (l + t)- 1 . 1 (3.10) 

Thus we proved (I3.3p 2 . 
Now we compute (13. 5J) 3 : 

• Interaction of transonic boundary layer solution and 3-rarefaction wave: 

/ I (V r:i -v*)\ + \ Vp (V b - «.) I di 
Jr + 

/oo 
V b (v*-V r *) +V? 3 (V b -v*) d£ 
l -a-+\ 3 (v*,e*)](l+t) 

= 0(l)5 b (l + 5 b t)- 1 

= 0(l)min{5,(l + £)-!} 

= 0(l)«jB(l + t)-5. (3.11) 

where in the first equality we have used the fact iv) in Lemma 2.3. 
Then we verify (13. 3j) ^ : 

• Interaction of 1-rarefaction wave and viscous 2-contact wave: 



14 



First we have 



/ \Vf(v m -V^)\d^ 



,_ + *l(Wm)| (1+( ) +QO 



. MH^gg) l(i + t) 



:— J5 + J6- 
Then we can compute 



j 5 = r ' \v^(v m -v r ")\dc 



r \-cr-+ Xliv ™' em) \{l+i) / n(t^„4\i 







1 + t 

C(l)5 d e- Ct , (3.12) 



and 

/>oo 

J 6 = / l^(«m- V ri )\d£ 

= 0(1)5, ' sup («m-V ri (*,0)=O(l)5« i e- ct . (3.13) 

g>[- --+ Al( "™ A " ) ](i+t) 

Similarly, we can estimate the interaction term 

/ |Vf (V d - v m ) I d£ = 0(l)5 d e~ ct . (3.14) 

</R + 

So (|3.3p ^ is verified. 
For (J33D 5 , that is 

• Interaction of 3-rarefaction wave and viscous 2-contact wave, which can be done 
similarly as (|3.3|) a. we omit the details for simplicity. 

Finally, we prove (I3.3p fi : 

• Interaction of 1-rarefaction wave and 3-rarefaction wave: 

Since Vp > 0, Vp < and the facts iii) and iv) in Lemma 2.3, one has 

/ I Vp {V r * -OI + I V[ 3 (y 1 - v m ) I d£ 

</R+ 

/+00 
vp (v* - V ra ) d£ 
i -a-+X 3 (v*,e*)](l+t) 

= 0(l)5 ri e- ct = 0(l)Se- ct . (3.15) 

Thus we justified (13.31) . The proof of (13. 4p can be done similarly, but the decay rates 
with respect to the time t may be higher. Therefore, we complete the proof of the wave 
interaction estimates in Lemma 3.1. ■ 

With the wave interaction estimation Lemma 3.1 in hand, we have the following 
Lemma: 
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Lemma 3.2. 

\\G(t)\\ L i + \\H{t)\\ L i =0(1)^(1 + 

\\G(t)\\ + \\H(t)\\= 0(1)5(1 + t)-\ (3.16) 

Proof. We can compute 

G l = | (P - P b - P ri - P d - P ra ) | 

= 0(l)\V b \ (\V r > - v *\ + \V d - v m \ + \V r > - v*\) 
+0(l)\V d \(\V b - v*\ + \V ri -v m \ + \V rs -v*\) 
+0(l)\vp\(\V b -v*\ + \V d -v m \ + \V r > -v*\) 
+0(l)\Vp\(\V b -vJ\ + \V^ -v m \+ \V d -v*\). (3.17) 

Thus by the wave interaction estimation Lemma 3.1, we have 

HdlUi =0(1)<J5(1 +*)-!§. 

1 13 

Similarly, ||ifi||x,i = 0(1)5» (1 + t)™ can be obtained. 

Now we estimate ||G 2 ||li and H^IU 1 - Note that in G 2 , besides the wave interaction 
terms, there are the error terms due to the i— rarefaction waves (i — 1, 3). So we can write 
G 2 as 

v i=l,3 7 ? v i=l,3 7 ? 

:= G21 + (j22- 

Since the wave interaction terms G21 can be verified similarly as G±, we only compute the 
error terms G22 due to rarefaction waves. 

\\G 2 ^=o(i)Y,(\\m\^ + \m^ v P)\\ 2 ) 

i=l,3 

= 0(1)55(1+*)-^ 

if we choose g > 14 in Lemma 2.3. 

In if 2 , besides the wave interaction terms and the error terms due to the %— rarefaction 
waves (i = 1,3), there exists the error terms H d due to the viscous 2— contact wave. We 
can compute that 

\\H d \\ L i =0(l)6 d (l + 0- 2 / r exp(- g ' ( * + *- t)2 )de 
= 0(1)5(1 + *)"§. 

The estimation of ||G|| and \\H\\ can be done similarly, thus the details are omitted. ■ 
3.2 Reformulation of the Problem 

Put the perturbation ((f), ip, £) around the superposition wave (V, [/, ©)(i, by 

(0, = (u, «, o)(t, - (f, c/, e)(t, 0, (3.i8) 
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then by (jl.7p and (13. ip . the system for the perturbation (<f>, if), i?)(t, £) becomes 



<p t - a_0 5 - ^ = 0, 

iP t -a^ + (p-P)^^ ~G, t>0, £>0, 

- + fa - PU t ) = k (* - + - <£f\ - H, 

(ip ,ipoM(0 ■■= (&^M)0U) (0,0,0), as e^+oo, 

(0, £ = o) = (n f/ rf , e rf )(t, e = o) - (v m , u m , e m ). 



Define the solution space X(0, T) to the above system by 



(3.19) 



X(0,T):=<j (<f>,if),$)(t,£) ((f ) ,iP^)eC([0,T};H 1 ), 5 G L 2 (0, T; L 2 ) , 

&L 2 (Q,T;H l ), N(T)=: sup ||(0, V, 0)(t)||i < e }, (3-20) 



Hereeo<| mm< l inf V"(t, £), inf 6 (£, £) > is a suitably small and positive constant 

I R.+ x R,-|- R.+ x R_|_ I 

to be determined. 

Since the proof for the local existence of the solution to the system (I3.19P is standard, 
the details are omitted. To prove Theorem 2.1, it is sufficient to prove the following 
a priori estimate by combining the local existence of the solution and the continuation 
process. 

Proposition 3.1 (A priori estimate) Let (<p,if>, i?) G X(0,T) be a solution to the 
system (\3.19§ in the time interval [0, T) with suitably small Eq, and the conditions in 
Theorem 2. 1 hold. Then there exist a positive constant C independent of T such that 



IK0,tM)(f)||?+ / [ll^(-r)|| 2 + ||(^,^)(r)||?] dr 



| J(U},U?,U?){i/>, 0)(t) \\ 2 dr <c(\\ (0o, V>o, #o) + ^ 



(3.21) 



3.3 Energy estimates 

To prove Proposition 3.1, we need the following several lemmas. First we give the following 
boundary estimates whose proof can be found in [21J. 



Lemma 3.3 (Boundary Estimates) [21J There exists the positive constant C such that 
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for any t > 0, 



|(0,<M)(T,O)| 2 dT<C5, 

( | We | + | W € | ) (r, 0) rfr < + /* || C ) (r) || 2 dr. 

./o 

(|0 T V|+0l)(r,O)dr<C5 + e / ||^(r)|| 2 rfr + C e / U^fdr, 

Jo Jo 

{\A^\+^!)(r,0)dr<C8 + e f \\^(r)\\ 2 dr + C t [ \\^(r)\\ 2 dr, 

Jo Jo 

(|0 T ^|+#?)(r,O)dT<Ctf + e / ||%(r)|| 2 dr + C £ / ||^(r)|| 2 rfr, 



'0 JO JO 

where e > is a constant to be determined and C e is the constant depending on e. 



Lemma 3.4 Let (0, 1/),$) G X(0,T) be a solution to the system ( Iff, iff) for some positive 
T and suitably small Eo > 0, and the conditions in Theorem 2.1 hold. Then there exist a 
positive constant C such that 



,^mm+ / \\MrW+\\(^mr)\\idr 



+ 



\y/(U>,U?,U?){<l>,0)(T)\\ 2 dT 



< C 



0,^,00)11? + ** / {l + r)-n{^Mr)\\ 2 dr 



\{<P,d)\ 2 d^dr. 



Proof. Step 1. Define 

$(77) := 77 — In 77 — 1 . 
Under the a priori assumption, there exist a positive constant C such that 

cry < $(77) < Crj 2 . 

Let 



(3.22) 

(3.23) 
(3.24) 



, v V 

Then a complicated but direct computation gives 



E t -F, + ^ + — ^ + P(C/| + + Up) 



v8 2 



V 10 



(3.25) 



Q, (3.26) 



where 



-PUf 



TV"®! 
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+ 
+ 

i=4 



nUzWt 2fiU^ kSS^ G{&&£ t*(Uz) 2 <lr& k(Q{) 2 

H 7, 1" T r nn \~ K 



vV 

e 



v9 



vve 2 



H 



vve vve 2 



v6 2 
V, 



=:^g,. (3.27) 

i=i 

Integrating (13. 26j) over [0,t] x R + yields 

II (0, ^ *) II 2 + f II ^)(r) frfr + jf || y/(Ul,U?,U?)fa *) 
< C||(0o,^o^o)|| 2 + C / |F(r,£ = 0)1^ + ^4 



fr)|| 2 rfr 



(3.28) 



8=1 



where = 0(1) / / Qid^dr. 
Jo Jn + 

From the boundary estimates in Lemma 3.3, we have 

' IFCr.e = 0)|dr < C5 + C5 f \\(^)(r)\\ldr. 
'o Jo 



We can compute that 



\{<f>^)\ 2 d^dr 



(3.29) 



(3.30) 



and 



h < C 



$)(T)\\ L ~(\\G(T)\\ L , + \\H(r)\\ Ll )dT 



<C6* / (l + r)-i6||(^,^)(r)|^||^^)(r)||^r 



< 6 jf || (^, € ) (r) frfr + C e 8 l e (l + J\l+r)- 



tf)(r)|| 2 dr (3.31) 



where and in the sequel e > is a small constant to be determined and C t is the positive 
constant depending on e. 

Now we calculate I3. By Cauchy inequality, we have 



t rt 
2 



h <e / 11(^,^)11^ + C e / / \(U^)\-UM 2 didT. 

'0 </0 ./R+ 



By Lemma 2.1-Lemma2.3, one has 

|(^,e 5 )| 2 <c 

By the techniques in [19 



+ 1 + + 5(1 + t)- 1 exp ( Gd{i + ^ 



(i + 50 4 



1 + t 



(3.32) 



(3.33) 



i/(«,oi = i/(*.e=o)+ / mm\ 

Jo 

< \ f{t,Z = o)\ + y/z\m\, 
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we have 




-\{^d)\ 2 didr 



R + (1 + *0" 



< C5 3 / |(0,tf)(r,£ = O)| 2 dr + C 



IK^.^II 1 



5^ 



dr 



<C6(1+ ||(^,^)(r)|| 2 dr 



Substituting ([3T33]) and (l3T34j) into (1332]) yields 



(3.34) 



h<C(e + 5) \\{^){r)\Ydr + C5 \\Mr)\\ 2 dr 



+ C6 + C6* J (l + r)-5||(0,^)(r)|| 2 rfr 



\{ct)^)\ 2 didT. (3.35) 



Then we have 



h = o(i) f f |(e^y/,c/|,e|,iy)||(0^)| 2 ^T. (3.36) 

Jo Jn + 

So I4 can be estimated similarly as I2 and 1%. 

Combining 1KT25]) . (13301) . 1EE5I1) . (I3~32l) . (13351) and (13361) . and then choosing 5 and e 
suitably small yield that 



||(0,V,tf)(t)H 2 + / ||(^,^)(r)H 2 dr+ / lly^^fX^X^II 2 ^ 



<C Ml (0o,^o^o) lr + f +C5b / ||(^,V«,%)(r)|| 2 dr 



+ / (l + r ) 



^,^)(r)|| 2 dr 



+ C5 




(l + r)" 1 esq) ( _ ^ + g - r )^ ) |(^)fdfr/-. 



JR 4 



l + r 



(3.37) 



Step 2. Differentiating (13.191b w.r.t. £ and multiplying it by |f yield 



2v 2 



0| \ _ 

2v 2 J £ v 3 v 2 



(3.38) 



Multiplying (13391) ■> by ^ gives 



(3.39) 
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H x (T3T38|) - f l339|) gives 



a. 



V 



V 



v 



^r^-^-^^ + Z^ 3 t7 

f v 6 v 6 \ vV J ^ v 



(3.40) 



Integrating ( I3.40P over [0, t] x R + , using the boundary estimations in Lemma3.3 and 
choosing 5 suitably small yield 



ll<M*)l| 2 + WMrWdr 



<C ||(V>o,0o f )|| 2 + 56 )+C5s / (l + r)^||(0^^)(r)|| 2 rfr 



C [5-«+e \\(^^)(r)\\ 2 + CMdr)\\ 2 \dr 



+ C5 f [ (l + r)- 1 expf- C, ^ + (7 - r)2 ) |(0^)| 2 ^ 

JO JR + V 



1 + T 



(3.41) 



Step 3. Multiplying f l3TT9l^ by then 



(p - + 



^(3.42) 



Integrating (I3.42p over [0, t] x R + yields 



+ / \\Mr)\\ 2 dr 



<C || (0o,^o^o) ||? + 5« )++C5-s / (l + r)-^||(0^^)(r)irrfr 



+ J {c(^+ e )||(^,%)(r)|| 2 + a||^(r)ir}rfr 

+ C6 Io L {1 + TrleXP {~ Cdi \ + +r T)2 ) (^ + d£dr 
where we use the following estimate 

ft \<P^M d idr< r||^(r)||||^(r)||||^(r)||^rfr 

Jo Jk+ Jo 



(3.43) 



< / ll^(r)||||^«(r)||'||^(r)||'dr 
Jo 

<e fu^r)\\ 2 dr + C e el f U^fdr. 
Jo Jo 



(3.44) 



Multiplying fET[9]U by then 
1? 



7-1 



0| 



0"_ 



2 ^ 



■0 



K 



+ -d 
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Integrating (I3.45P over [0,t] x R + yields 



K) 2 (u,y 



V 



+ H 



r). 



(3.45) 



%(r)frfr 



<c(||(0o^o^o)|| 2 + ^)++^ / (l + r)-i||(0^^)(r)|| 2 rfr 



+ 
+ C5 



o 



C tfi + c ||(^,%)(T)|| 2 + a||^(r)r 



'1 + r) exp 



'o JR+ 

where we use the following estimate 



C d {tj + <r-T) 
1 + r 



|(0,tf)| 2 d£dr, 



(3.46) 



JR 



<e ||(^,%)(r)||^r + C e ^ / ||(^)(r)||^r. 



(3.47) 



Combining (13.371) . (I3.4ip . (13.431) and (13.461) and choosing 5, e and e suitably small, we can 
complete the proof of Lemma 3.4. ■ 

Now to close the a priori estimates, the remaining thing is to compute the last term 
in the right-hand side of f!3.22j) which comes from the viscous contact wave. Here we use 
the method of the heat kernel estimation invented in [2]. 

Lemma 3.5. |2] Suppose that h{t,^) satisfies 

/i G L°° (0, T; L 2 (R + )) , ti£ G L 2 (0, T; L 2 (R + )) , h t - a_h 6 G L 2 (0, T; #-^+$,48) 
then 

2a(£ + ( x_r) 2 " 



1 + r) *exp 



o Jn 



1 + r 



/i" al^alr 



<C a ||/,(0,-)|| 2 + / /i 2 (r,0) + ||^(r,-)|| 2 + <^-a_^,( W a ) 2 / i > H _ lxH1 drfe.49) 



where 



w a (t,0 = -(l+tr* 



exp 



1 + t 



(3.50) 



and a > is a constant to be determined. 

Based on Lemma 3.5, we have the desired estimates in the following Lemma. 



Lemma 3.6 There exist a uniform constant C > such that if 8 and Eq are small enough, 
then we have 

jTjf (l + r)- 1 exp(- ^ i + a r - T)2 ^ |(<^,d)| 2 d£dr 
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<C ||(0o,Vo^o)||? + P )+CS-s / (l + r)-^!!^^^)^)!! 2 ^. (3.51) 



Proof. Step 1. First, let 



h = Pcj> + 



R 



7-1 



-x) 



(3.52) 



in Lemma 3.4. Then we only need to control the last term of (I3.49P on the right hand 
side. 

We have from the energy equation (I3.19p 3 . 



ht - <r_/i € = (P - p)^ + U^P -p)+ (P t - a.P^ 



(3.53) 



Thus 



/ (K- (J-h^(w a ) 2 h) H _ lxH1 dT 
Jo 



K 




JRj 



+ 




JR 4 



[(P - p)^ + U&P -p) + {P t - a_P ? )0 
_ &£.)-H] {w a fhdidr. 



(3.54) 



Notice that 
IK(t)IU<» < C a , = (1 + ^exp 



eg + <r_t)' 
1 + t 



-i 



(3.55) 



thus to control terms on the right hand side of (I3.54p . we only consider the term (w a ) 2 (P — 
p)hip£. By using the mass equation (I3.19p i and the momentum equation (13.19p 2 again, 
we have 



-fP(w a ) 2 h 
2v 



(7- 1)K) 2 ^ 2 



jP(w a ) 2 h(p 2 - 2( 7 - l)(w a ) 2 0/i 2 
2^ 

,a\2l, j,2 



7 P(w a ) 2 /i0 2 - 2(7 - l)(w a ) 2 0/i 2 



2?> 



7 P/i0 2 - 2(7 - l)0/i 2 



-fP(w a ) 2 h(p 2 -2h -l)(w a ) 2 (ph 2 , , rrN 

+ 2 ; 2 — + 



U ^ (Pt-a-Pz) 
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(w a)2 [4 ( 7 _l)/,_ 7 p0]0 

+ — [ht-a-ht). 



(3.56) 



Now the terms in the right hand side of (I3.56P can be estimated directly and in particular, 
we have 




o -/R 4 



< C 




o Jn^ 

t 



7 P(w a ) 2 /i0 2 - 2(7 - l){w a ) 2 (ph 2 
2V 2 

3 + |tf| 3 ) d£dT 



<c ||(0^)|ILII^IIII(0^)II^ 



<Ce 2 11(^,^,^)^)11^. 



(3.57) 



The other terms can be controlled by the similar procedure as Step 1 of Lemma 3.4. Thus 
the combination of the above estimates and Lemma 3.5 yield 

<cJ\\{Mo,#^\ + 8*)+C a 8* [ (l + r)-i||(0,^)(r)|| 2 dT 



+^)/7( 1+ ^ P (-^) 



|(0,tf)| 2 d£dr. (3.5* 



Step 2. Let 



cxp 



Ay 2 



l + t 



dy, 



(3.59) 



where A > is a constant to be determined 
Then 

From the fact p- P = m ~ p ^ , we have 

(m - P0) c vt(m - P(f>) 



\W t A - (T-Wf \ < C A (1 + i)'*. (3.60) 



- (ip t - a_^ c ) + y, 



u i _U i 
v V 



G. (3.61) 



Multiplying f |3~Ml) by W A (Rd - P<j>) implies 

w A {m-P(t)f\ w A {m-P4>f w\(m-P(j)) 2 



2v 



2v 



-W J 



2v 2 

{Rd-P(f>). 



(3.62) 
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Note that 



W A (i) t - a_il>z)(M - P<j>) = {W A i){R& - P<p)} t - o-{W A i)(M - P(j))} i 

-ij)(m - P<f>) (W t A - a^W^) 
-W A ip{(m - P<j)) t - a_{M - P0) ? }, 



(m - P(f>) t - <t-{r& - p<t>)t 



(7-i)U^-pK + Mv-f )+\ r 



K) 2 (Utf 



V 



H 



--yPi(> ( -(P t -(T-Pt)<f> 



(3.63) 



and 



we have 



where 



(3.64) 



■{(R&- P(j)) 2 + 7t»P^ 2 } = -{W A ip(R$ - P(f))} t - + Q A , (3.65) 



A . _ w\m - + lpwA 2 _ A _ (u, u 4 



E A : 



2v 

( 6 

aJW^iRd - P(j>) - (7 - 1)kW a i/) — - tt 

\ v V 



(3.66) 



and 

QA _ W\(P-p) 2 



+ (w A - a_w^) [m - p<p)ip - w A G{m - p<f>) 



+ W A i(){ (7-I) 



(P t -a-Pt)<l> + 



jp^A 
2 J 



- ,{W A (W - P0)} ? - I) - ( 7 - 1)k(W a ^ - f ) • (3.67) 



First, we have 



/ P A (r,0)c/r 



<C A 5 + C A 5 / ||(^,^)(r)||?dr. 
'o 



(3.68) 



The estimations of the terms concerned with W A are similar to those in Step 1 while the 
other terms are similar to those of Step 1 in the proof of Lemma 3.4. Thus integrating 
fl3"65D over [0,t] x R+ yields 



[ jf (1 + r)- 1 exp (- ^^l ) {(m - P0) 2 + ^} d^dr 



< Cj 



o, ^o, 0o) II? + ^ J + ^ / (! + ^ II (6 ^' 
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+ C A (6 + e ) jf (l + r^exp (- g^l±^£ ) tf)| 2 d£d<r. (3.69) 

Step 3. Combining ( 13. 58ft and ( 13. 69 p . then choosing A = 2a = Cd and setting 5, Eq 
suitably small, we can complete the proof of Lemma 3.6. ■ 

Proof of Proposition 3.1. Choosing 5, Eq suitably small in Lemmas 3.4 and Lemma 3.6, 
then using Gronwall inequality yield Proposition 3.1. ■ 
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